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ABSTRACT 


This  report  is  an  expansion  of  SAAC  Report  No.  4. 
The  SAAC  long  period  processing  methods,  based  on 
frequency-wavenumber  analysis,  are  capable  of  analyzing 
all  the  events  reported  in  the  LASA  and  NORSAR  Daily 
Event  Summaries  within  real-time.  The  capability  of 
LASA  to  detect  Rayleigh  waves  is  azimuth  dependent.  A 
comparison  of  NORSAR  and  LASA  on  a  data  set  of  Kurile 
islands  events  resulted  in  a  detection  of  surface  waves 
from  all  events  with  m^  4.0  as  reported  by  the  LASA 
Daily  Event  Summary.  A  study  of  central  Asia-western 
Soviet  Union  events  shows  that  NORSAR  M  values  are 
consistently  higher  than  those  at  ALPA  or  LASA,  the 
lowest  level  detections  being  made  were  on  events  with 

Ms  =  2*5  and  2*7  for  the  western  Soviet  Union  and 
central  Asia  respectively. 
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I .  INTRODUCTION 


This  report  describes  the  long  period  processing 
methods  developed  and  applied  in  the  SAAC.  Using  these 
methods,,  a  preliminary  evaluation  is  given  of  the  capa¬ 
bility  of  the  three  large  LP  arrays,  LASA,  NORSAR  and 
ALPA,  to  detect  and  measure  the  amplitude  of  Rayleigh 
waves  associated  with  P-wave  detections  reported  in  the 
LASA  and  NORSAR  Daily  Event  Summaries.  The  phrase 
"preliminary  evaluation"  is  deliberately  used  because 
long  period  data  are  now  rou'-.inely  analyzed  from  all 
three  arrays  and  so  the  data  on  which  our  evaluation 
is  based  are  continually  increasing.  Also  our  process¬ 
ing  methods  are  continuously  revised  and  improved  with 
a  resulting  improvement  in  the  capability  to  detect 
small  signals  and  separate  interfering  waves.  Thus  the 
results  and  conclusions  presented  in  this  report  must 
be  regarded  as  a  current  evaluation  which  will  probably 
be  modified  as  more  detection  results  become  available. 
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II.  LP  DATA  PROCESSING 


The  choice  of  data  processing  methods  was  influenced 
by  the  fact  that  a  large  number  of  events  have  to  be 
processed  in  real  time.  Time-domain  beamforming  and 
plotting  is  a  slow  process  which  suffers  the  disadvant¬ 
age  that  the  propagation  azimuth  and  velocity  must  be 
assumed.  Any  disturbance  visible  on  a  particular  beam 
may  be  associated  with  adjacent  beams  or  a  side  lobe. 

For  these  and  other  reasons,  a  frequency-wavenumber 
(f-k)  approach  has  been  developed  for  routine  analysis 
of  the  LP  data.  A  description  of  the  f-k  analysis  method, 
as  used  in  the  SAAC,  is  given  in  Appendix  I.  A  complete 
description  of  the  f-k  processing  system  is  currently 
being  written  as  SAAC  Report  No.  9. 

A  flow  diagram  of  the  process  is  shown  in  Figure  1. 
All  the  LP  data  from  NORSAR,  ALPA  and  LASA  are  recorded 
simultaneously  on  the  SAAC  low— rate  tape  each  tape 
containing  approximately  12  hours  of  recorded  data  for 
each  array.  The  data  are  used  directly  from  the  low- 
rate  tape;  no  intermediate  or  event  tape  is  used.  A 
particular  time  window  is  selected,  based  on  the  if-wave 
arrival  time  at  LASA.  The  individual  channels  are 
degl itched  and  edited  with  noisy  and  dead  channels 
discarded  automatically.  The  signal-to-noise  ratio 
detection  threshold  can  be  set  at  any  level  and  the 
velocity  range  of  interest  can  also  be  set  arbitrarily.  . 
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LIST  DETECTIONS  AS  FUNCTION 
OF  AZIMUTH.  VELOCITY.  PERIOD 
SIGNAL  POWER  AND  S/N 


Figure  1.  Schematic  illustration  of  the 


SAAC  LP  processing  program. 


The  multichannel  data  are  then  Fourier  transformed. 
Maxima  of  power  in  three-dimensional  f-k  space  are 
automatically  picked,  and  if  they  exceed  the  detection 
threshold  and  are  in  the  specified  velocity  range, 
these  maxima  are  listed  with  the  corresponding  azimuth, 
velocity,  period,  signal  power  estimate  and  signal-to- 
noise  ratio  estimate.  An  option  is  then  available  to 
detect  smaller  signals,  in  which  the  maximum  power  peaks 
and  associated  side  lobes  are  removed  and  second  maxima 
are  picked.  If  any  of  these  second  maxima  satisfy  the 
detection  criteria  they  are  also  listed  as  functions  of 
azimuth,  etc.  The  time  window  is  then  moved  forward 
or  backward  in  time  by  any  specified  amount  and  the 
whole  process,  including  the  editing,  is  repeated. 

The  decision  whether  a  peak  corresponds  to  the 
expected  Rayleigh  wave  is  based  on  several  criteria. 
First,  and  most  importantly,  the  time  must  be  correct. 
Using  a  256-second  window  and  25%  overlap,  a  time 
resolution  of  approximately  two  minutes  is  obtained  for 
a  particular  period.  This  resolution  enables  the 
dispersion  to  be  observed  and  this  itself  is  a  basis  for 
decision.  Second,  the  azimuth  and  phase  velocity  must 
be  correct.  Both  of  these  quantities  are  dependent  on 
period  so  the  expected  values  for  a  particular  region 
and  a  particular  array  are  obtained  by  observing  large 
amplitude  signals.  For  low  signal-to-noise  ratios  the 
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..oise  can  displace  the  location  of  a  peat  i„  f-k  space  ' 
slightly  and  hence  affect  the  velocity  and  azimuth  value 
listed.  For  such  cases,  deviations  from  the  expected 
azimuth  of  up  to  15°  are  tolerated.  Also,  the  analyzed 
data  are  examined  for  time  intervals  immediately  before 
and  after  the  detected  arrival  to  see  whether  the  peak 

is  persistent  or  spasmodic,  if  spasmodic-;  the  peak  is 
discarded  as  noise.  , 

1  :  t 

The  surface  wave  magnitude  Ms  is  derived  from  the  ' 
estimate  of  the  signal  power.  The  details  of  how 'this 
is  actually  done  are  described  in  Appendix  II .  Time- 
domain  beamforming  and  plotting  programs  are  available 
for  visual  display  of  the  recorded  data.' 

Programs  are  also  available  to  measure  this  spatial 
coherency  of  signals  within  the  area  of  an  array.  The 
spatial  coherency  gives  a  quantitative  Measure  of'  signal 
degeneracy  at  any  given  frequency  and  is  also  diagnostic, 
of  the  mechanism  causing  the  degeneracy  (Mack  and  Flinn, 
1972).  a  detailed  study  of  the  spatial' coherence  of 
Rayleigh  waves  at  the  three  large  arrays  is  currently 
being  written  and  will  appear  as  SAAC,  Report  No.  8 


I 


I 
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III.  CONTINUOUS  ANALYSIS  OF  LA£A  LP  DATA 

I 

In  SAAC  Report  No.  4,  an  experiment  was  described 
,  in  which  50  days  of  LASA  LP  recordings  were  analyzed. 

The  recording  period  was  in  the  summer  months  of  July 
and  August.  The  results  are  detailed  in  that  previous 

1  I1 

report  but  fbr  the  sake  of  completeness  the  main  con¬ 
clusions  will  be  outlined  here.  The  study  showed  that 
the  capability  of  LASA  to  detect  Rayleigh  waves,  using 
events  listed  in  the  paily  Event  Summary,  is  azimuth 
d^perident.  LASA  is  relatively  insensitive  to  events 

!  I 

in  South  America  compared  with  events  to  the  northwest. 

■  .  .  i  i  '  1  ! 

This  is  a  function  of  the  respective  travel  paths  and 
is  not  caused  by  az^Lmuthally  dependent  inoise. 

i  i 

The  background  noise  has  a  major  coherent  component. 

4  1  I  i 

This  appears  as  discrete  trains  of  Raylbigh  waves  prop¬ 
agating  across  the  array  from  all1  directions .  Because 
of  the  relatively, short  window  used  in  the  analysis,  the 

coherent  noise  does  not  appear  isotropic.  Rather  a 

.  '  ,  .  ’  .  i  ,  ' 

particular,  window  will  detect  one  of  two  of  the  discrete' 

v  I  ' 

wave  trains  andi  then,  as  the  window  is  moved  ,in  time, 

.  ,  i  -  i 

detections  will  be  made  from  other  directions.  The  50 

i  1 

days  of  analysis  showed  that  LASA  receives  a  good  deal 
of  coherent  noise  from  azimuths,  around  60°  to  70°.  ,  A 
good  deal  of  eneirgy  lies  in  the  microseismic  band  within 

t  I  .  I 

the  period  range  of  16  to  18  seconds.  However,  some 

!  I 

appear  with  periods  of  around  20  seconds.  Hence  the 

I  ' 
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threshold  of  detections  for  Rayleigh  waves  is  quite 
dependent  upon  the  nature  and  direction  of  coherent  noise 
at  any  particular  instant.  The  microseisms  from  the 
above  mentioned  azimuths  are  considered  to  be  associated 
with  sea  wave  activity.  Triangulation  using  the  other 

arrays  usually  puts  the  main  source  in  the  Labrador 
Sea. 

Another  important  outcome  from  this  study  was  the 
effect  of  a  large  earthquake  on  the  detection  capability 
of  the  array.  An  earthquake  near  Honshu  with  =  6.8 
produced  such  large  surface  waves  that  no  further 
detections  were  made  for  approximately  four  hours. 

After  the  initial  Rayleigh  wave  had  passed,  scattered 
waves  were  detected  from  a  variety  of  azimuths  and, 
even  though  much  smaller  than  the  main  arrival,  they 
were  sufficiently  large  to  mask  Rayleigh  waves  from 
smaller  earthquakes.  Presumably  NORSAR  and  ALPA  were 
similarly  affected. 
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IV.  NETWORK  EVALUATION 


Kuriles  Study 

SAAC  Report  No.  4  gave  some  results  of  a  LASA-NORSAR 
comparison  using  events  in  the  Kurile  Islands  region, 
an  area  which  is  equidistant  from  the  two  arrays.  This 
study  is  continuing  and  the  current  state  of  the  eval¬ 
uation  is  presented  in  this  report.  Now  that  ALPA  is 
working  consistently,  events  from  the  Kuriles  recorded 
at  ALPA  are  also  being  routinely  processed  and  these 
results  are  included.  A  listing  of  all  the  events  used 
is  given  in  Appendix  III.  This  listing  also  gives  the 
result  of  processing  each  event  for  each  array. 

Ms:m.  plots  for  the  Kuriles  events  are  shown  in 
Figure  2  for  LAS A  and  Figure  3  for  NORSAR.  The  two 
distributions  are  very  similar  except  for  the  fact  that 
the  NORSAR  M  values  are  generally  higher  than  those 
measured  at  LASA.  The  difference  in  travel  paths  could 
account  for  this.  Rayleigh  waves  received  at  LASA 
travel  through  the  north  Pacific  and  the  Aleutian  arc, 
then  through  western  Canada,  altogether  a  rather  complex 
path.  The  path  from  the  Kuriles  to  NORSAR  may  be  con¬ 
sidered  a  bit  simpler  as  the  bulk  of  the  propagation 
occurs  through  the  northern  part  of  the  Asian  continent. 
Marshall  and  Basham  (1972)  suggest  a  correction  of  +0.13 
to  an  Mg  value  measured  on  a  wave  propagating  over  a 

mixed  path  relative  to  the  M  value  measured  on  a  wave 

s 
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traveling  in  a  continental  structure  only.  This  suggestion 
generally  agrees  with  the  data  presented  in  this  report. 

It  is  of  interest  to  note  that  each  of  these  two  M  : 

s 

distribution?  may  possibly  indicate  the  two  distinct 
groups  suggested  by  Lacoss  (1971). 

Figures  4  and  5  show  a  breakdown  of  the  Kurile 
Islands  events  in  terms  of  detection  and  non-c(etection. 

All  events  reported  by  the  SAAC  in  the  period  May  - 
December  1971  for  which  4.0  <11^55.0  have  been  used. 

The  hatched  areas  denote  events  missed  due  to  the  pres¬ 
ence  of  another  larger  surface  wave  in  the  time  period 
when  a  Rayleigh  wave  from  the  area  of  interest  was 
expected.  These  interfering  waves  arrived  on  azimuths 
close  to  the  back  azimuth  from  the  Kurile  Islands  and  so 
array  processing,  even  the  "stripping"  mentioned  in  the  data 
processing  section,  was  not  of  any  help.  In  the  mag¬ 
nitude  range  4. 5  _<  m^jc  4. 9  the  three  events  missed  at 
each  array  because  of  interference  were  the  same.  In 
the  range  4.0<  m^  4.4  the  four  missed  at  NORSAR  were 
also  hidden  by  interference  at  LASA.  The  fifth  event 
hidden  at  LASA  was  detected  at  NORSAR.  Very  few  events 
were  missed  when  there  was  no  apparent  interference. 

For  body  wave  magnitudes  greater  than  or  equal  to  4.0, 
as  reported  in  the  SAAC  Daily  Event  Summary,  only  one 
was  not  detected  at  NORSAR.  This  had  =  4.1  and  it 
was,  in  fact,  detected  at  LASA.  Of  the  three  not 
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detected  at  LASA  two  were  detected  at  NORSAR  and  for  the 
other,  the  NORSAR  data  were  bad  on  the  low-rate  tape. 

If  we  accept  that  a  detection  at  just  one  cf  the  two 
arrays  is,  in  fact,  a  true  detection  then  apart  from 
one  miss  at  LASA  when  the  NORSAR  data  were  bad,  all 
events  with  "^>4.0,  as  reported  by  the  SAAC  Daily  Event 
Summary,  were  detected  when  there  was  no  obvious  inter¬ 
ference.  it  would  appear,  therefore  that  for  SAAC  n^- 
4.0  tne  most  serious  cause  of  non-detection  of  events  in 
the  Kurile  Islands  is  interference  from  other  events. 
Approximately  15%  of  this  data  set  could  not  be  detected 
at  either  array  due  to  interference. 

Processing  of  events  with  SAAC  m^  4.0  is  currently 
being  done  in  order  to  try  to  establish  some  working 
threshold  for  the  Kurile  Islands  region. 
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Central  Asia-western  Soviet  Union 

ALPA  and  NORSAR  are  conveniently  situated  for 
studying  surface  waves  emanating  from  the  Sino-Soviet'  ■ 
region.  The  ray  paths  from  a  particular  event  to  each 
array  are  approximately  orthogonal  but  are  generally 
not  equidistant  except  for  a  region  covering  Lake  Baikal, 
eastern  Mongolia  and  the  eastern  part  of  China.  However, 
as  the  central  Asian-western  Soviet  Union  region  is  of 
particular  interest,  an  ALPA-NORSAR  comparison  has  been 
made  using  events  from  that  area  rather  than  from  a  more 

localized  source  region  such  as>  the  Kuriles  an  the  LASA-  < 
NORSAR  study. 

The  data  from  ALPA  has  been  generally  good  only  since 
October  1971.  The  NORSAR  data,  transmitted  via  the  TAL, 
has  improved  towards  the  end  of  that  year  also.  The 
study  has,  therefore,  been  essentially  confined  to  events 
occurring  fairly  recently.  Another,  and  more  serious 
qualification  on  this  evaluation  is  that  much  of  the 
region  in  question  is  at  or  beyond  the  teleseisfnic  limit 
of  P-wave  detection  for  LASA.  Thus  LASA  has  a  relatively 
high  threshold  of  detection  for  events  in  that  area  and 
the  body  wave  magnitude  measurements  can  differ  quite 
drastically  from  those  reported  by  NOAA.  The  NORSAR 
Daily  Summary  is  not  edited  in  a  similar  manner  to  the 
LASA  Daily  Event  Summary,  so  we  have  to  refer  to  the 
NORSAR  bulletin  which  is  the  NORSAR  EP  output  trans- 
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mitted  via  the  TAL;  this  has1  been  available  only  recently. 
Presumably  this  bulletin  has  the  same  shortcomings  as  the 
SAAC  EP  output  bulletin  (SAAC  Report  No.  1)  and  hence 

I  '  :  ! 

must  be  used  with  caution.  However,  for  events  in  the 
central  Asian  region  the  body  wave  magnitudes  reported 
by  NORSAR  tend  to  be  more  in  agreement  with  the  NOAA 

values  and  so,  wherever  possible,  these  values  are  taken 

\ 

in  preference  to  the  LASA-SAAC  values. 

1  :  I 

The  list  of  events  used  is  tabulated  in  Appendix  III 
together  with  the  Mg  values  at  each  array  if  available. 
More  measurements  were  obtained  from  NORSAR  because  the 
data  were  consistently  better.  Currently  this  situation 
has  improved  and  good  data  are  being  received  from  both 

i 

arrays  at  comparable  rates. 

Although  these  events  are  by  no  means  a  complete 
data  set,  no  particular  selection  rules  have  been  applied 
and  so  the  results  are  most  probably  representative  of 
a  larger  set.  Only  events  reported  by  NOAA  as  being 
deep  (  >  60  km)  have!  been  omitted.  Some  deep  events  may 
be  included.  Very  few  small  magnitude  events  (m^^.4.0) 
are  included  simply  because  they  do  not  pass  the  LASA- 
SAAC  EP  threshold.  This  situation  is  being  remedied  as 

i  1 

we. become  more  familiar  with  the  NORSAR  bulletin. 

Figure  6  shows  Mg:mjj  values  for  NORSAR.  Good  sep¬ 
aration  between  earthquakes  and  presumed  explosions  is 
observed.  Figure  7  gives  the  equivalent  values  for 


i  i 
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ALPA.  When  the  two  sets  are  directly  compared,  as  shown 
in  Figure  8,  a  significant  result  emerges.  In  all  cases 
for  earthquakes  the  Mg  value  calculated  at  NORSAR  is 
higher  than  the  same  value  at  ALPA.  The  minimum  sep¬ 
aration  is  0.2  Mg  units  and  the  maximum  is  0.8  Mg  units. 
Again,  as  in  the  LASA-NORSAR  comparison,  the  travel  paths 
are  geologically  different.  The  path  to  NORSAR  traverses 
stable  continental  and  shield  structure  whereas  Rayleigh 
waves  recorded  at  ALPA  travel  from  a  continental  source 
skirt  the  Arctic  Ocean  basin  and  pass  into  the  tectonically 
complex  Alaskan  region.  Spatial  coherence  studies  have 
shown  that  surface  waves  generated  in  central  Asia  exhibit 
very  little  variation  across  the  NORSAR  array,  whereas 
there  is  quite  a  significant  amount  of  variation  across 
ALPA.  During  quiet  periods  the  RMS  background  noise 
level  is  approximately  the  same  for  the  two  arrays  for 
beams  focused  towards  central  Asia  (Texas  instruments, 

Inc.,  personal  communication).  This  implies  that  NORSAR 
is  a  better  long  period  array  than  ALPA  for  detection 
of  surface  waves  generated  in  the  central  Asian  region. 

The  few  M  values  calculated  at  LASA  for  events  in  this 
area  indicate  that  NORSAR  again  has  systematically  higher 

values . 

In  Figure  9  a  composite  set  of  M^n^  values  is  given 
for  the  central  Asia-western  Soviet  Union  region.  All 
the  earthquakes  are  recorded  by  ALPA  or  NORSAR.  In  cases 
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Figure  9.  Combined  NORSAR-ALPA  Ms : 
from  central  Asia  and  the  western 
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Soviet  Union. 
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where  both  arrays  recorded  the  event  the  higher  M 

s 

value  (always  NORSAR)  is  used.  The  same  applies  to  the 
presumed  explosions,  the  higher  value  being  used  also. 

A  few  values  for  presumed  explosions  recorded  at  LASA 
early  in  1971  have  been  included.  These  are  all  docu¬ 
mented  in  Appendix  III.  A  good  separation  is  mainLdined 
between  earthquakes  and  presumed  explosions. 

The  lack  of  small  amplitude  events  would  make  it 

somewhat  presumptious  to  talk  in  terms  of  a  threshold 

of  detection.  However,  using  the  frequency  beamforming 

techniques  described  in  Appendix  I  an  Mg  value  as  low  as 

2.5  has  been  measured  at  NORSAR  for  a  western  Soviet  Union 

source  and  an  M  =  2.6  at  ALPA,  and  for  the  central 
s 

Asian  region  the  lowest  values  sj  far  measured  are  2.7 
and  2.9  for  NORSAR  and  ALPA  respectively. 

Although  NORSAR  is  appearing  to  be  a  very  good  LP 

array  for  events  from  the  east,  there  is  at  least  one 

potentially  serious  handicap.  Coherent  noise,  at  Rayleigh 

wave  velocity,  has  been  observed  propagating  from  an 

azimuthal  range  50°  to  90°.  In  one  particular  case  this 

noise  interfered  with  the  reception  of  a  Rayleigh  wave 

generated  by  a  presumed  explosion  in  the  estern  Soviet 

Union.  This  example  was  shown  in  SAAC  Report  No.  4. 

The  event  was  detected  at  LASA  and  ALPA,  a  fact  which 

establishes  the  value  of  the  network,  but,  for  smaller 

sources,  NORSAR  is  the  only  array  capable  of  making  a 

* 

detection,  owing  to  its  proximity  to  western  USSR. 
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Conclusions  1  ' 

The  current  analysis  indicates  that  NORSAR  is  the 
best  of  the  three  large  arrays  for  detecting  and  measur¬ 
ing  Rayleigh  waves  from  central  Asia  and  the  western 
Soviet  Union.  This  is  not  just  due  to  the  fact  .that  it 
is  closer  but  also  because  signals  show  virtually  no 
degeneration  across  the  array.  However,  the  occurence 
of  the  coherent  noise,  observed  at  NORSAR  propagating 
from  the  east,  establishes  that  at  least  a  .second  array, 
in  this  case  ALPA,  is  necessary  to  avoid'  mis  identification. 

i 

Analysis  of  small  events  is  now  in  progress  to  estimate 
a  network  threshold  but  so  far  we  know  only  that  signals 
with  Mg  =  2.7  in  central  Asia  and  2.5  in  western  Soviet 
Union  have  been  identified. 

J 

The  NORSAR-LASA  study  of  Kuriles  events  shows  that 
for  mb>^4.0  interference  from  events  in  adjacent  seismic  1 
regions  was  the  main  cause  of  hon-detection*  In  the  data 
set  used,  which  was  derived  from  the  LASA/SAAC  Daily  Event 
Summary,  all  events  with  ^>4.0  were  detected  when  no 
interference  was  present,  with  a  detection  at  just  one 
array  (and  not  necessarily  the  other)  being  considered 
a  true  detection.  The  wave  amplitude*’  at  NORSAR  were 
generally  larger  than  the  equivalent  amplitudes^  at  LASA 
A  simpler  propagation  path  and  better  signal  coherence 
at  NORSAR  is  considered  to  explain  this  difference.'  1 


» 


V.  CURRENT  AND  FUTURE  WORK1 

1  ! 

1  '  1 

The  current  LP  processing  system  can  analyze  the 
LP  data  from  all  three  arrays,  looking  for  and  measuring 
Rayleigh  waves  associated  with  all  events  reported  by 
the  LASA  and  NORfJAR  event  summaries,  all  at  a  real-time  1 
rate.  However,  this  complete  analysis  requires  an  in-r 
ordinate  amount  of  360/44  time,  considering  the  general 
computer  requirements  of  SAAC,1  and  so  feasibility  tests 
are  being  done  to  see  how  much  LP  analysis  can  be  per¬ 
formed  pn  a  routine  basis  using  both  the  360/44  and  the 

i  1  1 

,  360/40B.  Currently  all  events  from  the  Sino-Soviet  and 

adjacent  i  regions  reported  by  the  LASA  or  NORSAR  event 

,  1  1 

summaries  are  routinely  analyzed,  and  these  data  are 
being  collected  for  presentation  in  a  later  report. 

1  .  '  1  1 

A  comparison  study  between  f-k  analysis  and  a  cuasi- 

matched  filtering  scheme  is  being  initiated.  The  filters 
are  derived  from  composite  group  velocity  curves  gener¬ 
ated  for  any  source-receiver  combination.  This  particular 
patched  filtering  approach  is  beirtg  used  because  it  is 
considered  to  be  the  only  one  capable  of  handling  many 
events  routinely.  1  !  1 


I 
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APPENDIX  I. 


Frequency-Wavenumber  Analysis 

It  was  mentioned  in  the  text  that  frequency-wave- 
number  (f-k)  analysis  is  nothing  more  than  beamforming 
in  the  frequency  domain.  The  method  takes  advantage  of 
the  fact  that  the  signal-to-noise  ratio  is  frequency 
dependent  and  so  the  beamforming  is  done  frequency  by 
frequency.  Also  by  staying  in  the  frequency  domain  a 
great  many  beams  can  be  examined  rapidly,  in  practice 
this  means  that  the  azimuth  and  velocity  of  a  signal 
need  not  hs  assumed;  one  merely  accepts  the  beam  with 
the  maximum  power.  The  position  of  this  maximum  in  the  wave- 
number  domain  defines  the  azimuth  and  velocity  of  the  signal. 
The  process  is  defined  by  the  expression  (Smart  and  Flinn, 1972) 

-?  N  (1) 

P(f,k)  =  N“  I^I^An(f)  exp  [2tt1  (<}>n  (f)  -  k  •  rn)  ]  |  2 

where  An  (f )  exp  [ (f)  ]  is  the  Fourier  transform  of  the 
n'th  channel  located  at  position  vector  r  .  This  ex¬ 
pression  represents  a  three-dimensional  ^ace  with 
frequency  being  one  dimension  and  the  wavevector  k  being 
the  other  two.  For  computation  k  is  resolved  into  a 
Cartesian  coordinate  system  with  ky  related  to  geographic 
north  for  a  particular  array  and  kx  related  to  geographic 
east.  For  a  particular  time  window  this  space  is  searched 
for  maxima,  the  position  of  each  maximum  defining  the 
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velocity  and  azimuth  of  wave  motion  propagating  coherently 
across  the  array.  The  power  of  a  maximum,  P(f,k),  is 
the  signal  power  estimate  at  frequency  f,  and  the  signal- 
to-noise  ratio  estimate  is  defined  as: 

S/N  =  t-  ^  m 

-  pcf.k)  (2) 

where  Av(f)  is  the  average  channel  power  at  frequency  f. 
Obviously  this  estimate  is  better  for  high  signal-to- 
noise  ratios.  For  a  particular  array  geometry,  expression 
(2)  has  a  minimum  value  above  which  it  is  certain  that 
a  detection  was  not  made  via  an  array  side-lobe.  Below 
this  particular  minimum  value  a  detection  may  or  may  not 
have  been  made  by  a  side- lobe.  Expression  (2)  is  also 
indicative  as  to  the  goodness  of  the  estimate  of  signal 

power  P (f , k) .  The  higher  the  S/N  value  the  better  the 
estimate. 

An  example  of  how  effective  a  frequency  domain 
beamforming  technique  can  be  is  illustrated  in  Figure 
10.  The  beam  trace  of  a  Rayleigh  wave  originating  in 
the  Kurile  Islands  and  recorded  at  LASA  was  added  to 
LASA  noise,  from  the  same  beam.  The  amplitude  of  the 
signal  was  successively  lowered  but  the  noise  was  held 
at  constant  amplitude  and  at  each  step  the  data  were 
subject  to  the  described  analysis.  The  top  trace  in 
Figure  10  shows  the  noise,  with  an  RMS  value  of  6  my 
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(relative  to  the  25  second  period  calibration).  The 
signal,  with  a  peak-to-peak  value  of  12  my,  has  been 
added  but  is  quite  invisible.  This  is  now  equivalent  to 
an  Ms  value  of  2.6  as  recorded  at  LASA.  The  lower  trace 
is  the  same  data  filtered  with  a  bandwidth  of  20  to  40 
seconds  period.  The  signal  is  now  visible  but  is  dis¬ 
torted  by  the  additive  noise  and  is  indistinguishable 
from  other  wave  groups  which  have  been  passed  by  the 
beam.  However,  the  described  method  of  f-k  analysis 
assigned  the  correct  azimuth  and  velocity  to  the  signal 
and  assigned  other  azimuths  and  velocities  to  the  other 
wave  groups.  This  experiment  was  repeated  three  times 
using  different  noise  and  in  each  case  the  results  were 
essentially  the  same.  One  may  argue  that  similar  results 
might  be  obtained  in  the  time  domain  by  very  narrow  band 
filtering  and  forming  a  great  number  of  closely  spaced 
beams.  This  may  be  true,  but  the  computation  time 
involved  would  be  very  much  longer. 

As  was  mentioned  in  the  main  text,  the  interference 
of  two  separate  Rayleigh  waves  within  the  area  of  the 
array  is  a  sufficiently  common  phenomenon  as  to  warrent 
serious  consideration  in  any  processing  scheme.  Further¬ 
more,  using  the  relatively  short  time  windows  in  the  f-k 
analysis,  the  background  noise  appears  to  consist  mainly 
of  discrete  Rayleigh  wave  trains  arriving  from  many 
azimuths.  These  wave  trains  interfere  with  small  expected 
arrivals  in  much  the  same  manner  as  two  larger  events. 
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If  two  Rayleigh  waves  are  of  comparative  size,  f-k  ' 
analysis  usually  detects  the  presence  of  both,  and 
assigns  the  appropriate  azimuths  and  velocities  as  the 
time  window  moves  through  the  event,  because  it  is  highly 
unusual  for  the  two  waves  to  exhibit  identical  time- 
frequency  histories  in  the  array  region.  However,  signal 
power  estimates  may  be  in  error  even  thoiigh  the  events 
were  correctly  identified.  Also,  if  one  event  is  really 
dominant  then  it  may  completely  mask  the  smaller  arrival 
For  such  cases  a  filtering  scheme  has  been  devised 
(Smart,  personal  communication)  whereby  the  f-k  space  maxima 
are  removed,  along  with  associated  side-lobes,  allowing  the 
detection  of  secondary  maxima.  This  process,  referred  to  as 
stripping,  is  accomplished  by  removing  the  signal 
estimate  from  each  channel  in  the  complex  frequency 
domain,  frequency  by  frequency.  Specifically.,  the 
transform  of  the  j'th  channel  after  signal  subtraction, 

Fj(f),  is  expressed  as  1 

i 

F'j(f)  =  (f)exp[-27iik.r^  ]  i  | 

I 

1 

1  N  I  * 

JT  exp[27rikTj] 
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i  An  example  of  this  stripping  process  was  given  in 
SAAC  Report  No.  4.  In  that  particular  case  two  signals 

f 

were  added  as  though  recorded  at  ALFA  from  two  distinct 
azimuths  with  a  20  db  difference  in  signal  size.  It  is 
much  more  satisfactory  to  illustrate  the  process  with 
real  situations;  and  two  such  examples  are  given  in  this 
report.  The  first  example  illustrates  the  separation 
of  two  unequal  but  azimuthally  distinct  signals  recorded 
by  ALPA.  The  wavenumber  plane;  for  a  period  of  18.3 

,  S 

seconds,  is  shown  in  Figure  11.  An  arrival  from  an 

.1  o  ' 

azimuth  of  65  is  dominant  and  there  is  only  a  suggestion 

.  'I 

of  a  second  s.ignal  12  db  down  from  an  azimuth  of  187°. 
i  This  same  plane  is  shown  .in  relief  in  Figure  12  which 
effectively  illustrates  the  signal  power  difference. 

Figures  13  and  14  show  the  wavenumber  plane  and  the 
relief  version  after  application  of  the  stripping  process. 
The  second  signal  is  now  quite  clear  and  is  in  fact 
dominant .  ! 

l 

j 

The  second  example  shows  how  a  small  signal  was 
detected  at  ALFA  in  the  presence  of  another  small  signal  and 

i  ,  •  m 

i 

noise.  The  wavenumber  plane  and  relief  version  are 
shown  in  Figures  15  and  16  respectively.  The  signal  of 
interest  is  masked  by  anbther  one  close  by.  After 
stripping.  Figures  17  and  18 1  show  the  second  smaller 

I 

signal  quite  clearly.  This  small  signal  occurred  at  the 
expected  arrival  time  at  ALPA  of  a  Rayleigh  wave  from 
a  small  presumed  explosion  (SAAC  =  4.9)  in  the 

i  .  ,  • 
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western  Soviet  Union.  The  signal  amplitude  is  calculated 

as  11  my  corresponding  to  a  surface  wave  magnitude,  M  , 

of  2.6.  For  comparison  the  time  domain  beam  plot  is  S' 

Shown  in  Figure  19  and  it  is  obvious  that  a  pick  could 

not  be  confidently  made  at  the  expected  arrival  time 

because  the  signal  is  indistinguishable  from  other  dis- 
turbances. 
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Figure  11.  Wavenumber  plane  at  a  period  of  18.3  sec  for  ALPA  on 
4  October  1971  for  the  time  interval  1037  48.0  -  1042  3.0  GMT. 

The  suspected  smaller  signal  is  12  db  down  from  the  main  dominant 
peak. 
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WAVENUMBER  PLANE  AT  18.3  SEC. 

ALPA  4  OCT.  1971  10  37  48.0  -  10  42  3.0  GMT 


SUSPECTED 
SMALL  SIGNAL 


In3ngure'll?Uef  illustrati°"  °*  the  wavenumber  plane  shown 
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Fi9l?fe  i3,.The  wavenumber  plane  shown  in  Figure  12  after  removal 
of  the  dominant  signal  by  the  stripping  process. 
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Figure  15.  Wavenumber  ,plane  for  ALPA  on  the  19  September  1971  for 
the  time  interval  1125  00.  -  1129  15.0  GMT.  The  period  is  26.9  sec. 
The  smaller  hidden , arrival  to  the  north  ii  merged  with  a  slightly 
larger  arrival  to  the  southwest. 
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IDDEN  SIGNAL 


WAVENUMBER  PLANE  AT  2M  SEC. 

ALPA  19  SEPT.  1971  11  25  0.0  -  11  29  15.0  GMT 


Figure  16. 


Relief  version  of  the  wavenumber  plane  in  Figure  15 


ft  I  - 


T5e  ®ame  fave"ufber  Plane  as  in  Figure  15  after  removal 
interfering  signal  by  the  stripping  process.  The  second 
signal  is  now  distinctly  resolved. 
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WAVENUMBER  PLANE  AT  26.9  SEC. 

ALPA  19  SEPT.  1971  11  25  0.0  -  11  29  15.0  GMT  (After  stripping  ) 


Figure  18.  Relief  version  of  the  wavenumber  plane  shown  in  Figure  17 
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M  Measurement 

s _ 

As  was  stated  in  the  main  text,  one  reason  a  fre¬ 
quency  domain  analysis  scheme  has  been  devised  is  because 
of  the  large  volume  of  array  data  to  be  processed  in 
real-time.  However,  for  discrimination  purposes  the 
ground  motion  amplitude  of  the  surface  waves  is  desired 
in  order  to  calculate  the  surface  wave  magnitude  M  . 

This  is  conventionally  done  by  picking  the  maximum 
excursion  on  a  time-domain  plot.  To  avoid  confusion  the 
Mg  values  given  in  this  report  are  what  would  be  obtained 

by  measuring  amplitude  from  a  beam  plot  and  applying  the 
formula 

Ms  =  lo2  a/T  +  log  A  +  1 0 12 


where  the  symbols  have  their  usual  meaning.  The  signal 
power  estimate  for  a  particular  source  region,  array  and 
time  window  length  is  empirically  related  to  the  ground 
motion  amplitude  of  surface  waves  obtained  from  time- 
domain  beams  with  good  signal- to-noise  ratios.  To  illus¬ 
trate  how  this  relation  is  obtained,  a  large  amplitude 
Rayleigh  wave  from  the  Kurile  Isles  was  added  to  beamed 
noise  of  NORSAR  and  LASA.  A  series  of  f-k  analyses  were 
run  with  successively  reduced  signal  amplitudes,  and  the 
square  root  of  the  signal  power  estimate  was  plotted  as 
a  function  of  maximum  signal  amplitude.  The  results  for 
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LASA  are  shown  in  Figure  20.  A  good  linear  relationship  is 
obtained  for  higher  amplitudes  with  a  small  deviation 
as  the  amplitude  approaches  zero  resulting  in  a  finite 
power  associated  with  the  noise.  The  slope  of  the  line 
is  the  empirical  relation  between  surface  wave  amplitude 
and  the  square  root  of  the  signal  power  estimated  from 
the  f-k  analysis  procedure.  For  any  particular  region 
this  relationship  is  obtained  as  a  function  of  distance 
from  the  array.  To  illustrate  the  reliability  of  the 
method  Figure  21  shows  the  comparison  between  Mg  values 
measured  from  beam  plots  and  from  the  f-k  power  for  a 
set  of  Kurile  events  recorded  at  NORSAR  and  LASA.  The 
points  are  seen  to  cluster  around  a  line  of  unit  slope. 

The  small  amount  of  scatter  is  understandable.  For  one 
reason,  the  period  which  has  the  dominant  power  according 
to  the  f-k  analysis,  very  rarely  agrees  exactly  with 
the  period  of  the  maximum  amplitude  measured  from  the 
seismogram.  Another  is  that  the  signal  power  estimate 
varies  slightly  depending  on  the  position  of  the  window 
relative  to  the  maximum  amplitude  on  the  seismogram. 

This  effect  decreases  with  increased  overlap  of  time 
windows.  However,  the  trade-off  for  this  improvement 
is  increased  processing  time  for  a  given  record  length. 
Practically  speaking,  a  25%  overlap  is  found  to  be 
sufficiently  accurate  for  events  further  than  40°  from 
the  array.  For  closer  distances,  western  Soviet  Union 
to  NORSAR  for  example,  a  50%  overlap  is  used. 
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All  in  all,  the  conventional  and  f-k  methods  of 
measuring  Ms  agree  very  closely.  The  advantages  of  f-k 
are  speed,  automatic  selection  of  the  period  with  high¬ 
est  signal-to-noise  ratio,  and  the  certainty  that  the 
power  measurement  used  is  truly  from  the  beam  of 
interest. 
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Table  of  events  analysed  and  used  for  this  report. 
The  body  wave  magnitude  is  the  value  reported  by  the 
SAAC  Daily  Event  Summary.  Under  eaoh  array  name  is  the 
surface  wave  magnitude  Ma  if  detected  and  measured. 

A  star  <*)  signifies  the  Rayleigh  wave  was  not  detected 
and  the  letters  BD  indicate  bad  data  on  the  low  rate 
tape  at  the  expected  Rayleigh  wave  arrival  time.  If  the 
LP  data  were  not  processed  up  to  the  time  of  this  report 
a  gap  has  been  left.  The  letter  I  denotes  that  strong 
interference  from  an  adjacent  seismic  region  prevented 
any  possible  detection  being  made. 
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